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' If the causality condition [the speed of sound always remains less than that of light in 

i vacuum, i. e., V < c = 1] is imposed on the spheres of homogeneous energy density, the 

, 'ratio of the specific heats', 7 < 2.59457, constraints the compaction parameter, u[= 

"* **> ■ (M/a), mass to size ratio in geometrized units] of the dynamically stable configurations 

i-S-H ' < 0.34056 [corresponding to a surface redshift (z a ) < 0.771]. Apparently, The maximum 

value of u obtained in this manner belongs to an absolute upper bound, and gives: (i) The 
maximum value for static neutron star masses as 5AMq, if we substitute the density at 
f^ j 1 the surface of the configuration equal to the average nuclear density, E = 2x 10 14 g cm -3 

[e.g. Nature, 259, 377 (1976)]. (ii) However, if the density of the static configuration is 
constrained to the value 1.072 X 10 14 g cm -3 , by imposing the empirical result that 
the minimum rotation period of the fastest rotating pulsar known to date, PSR 1937 + 
21, is 1.558 ms, the maximum mass value for static neutron stars exceed upto 7.4Af(j. 
These masses have important implications for the massive compact objects like Cyg 
r% \ X-l, Cyg XR-1, and LMC-X3 etc., which may not, necessarily, represent black holes. 

(iii) The minimum rotation periods for a static 1.442Mq neutron star to be 0.3041 ms. 

(iv) A suitable stable model of ultra-compact objects [u > (1/3)] which has important 
astrophysical significance. 
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1. Introduction 

Incompressible fluid spheres of uniform energy density E in General Relativity were 
first discussed by Schwarzschild. 1 The importance of this solution in General Rel- 
ativistic stellar structures is apparent, because it gives an absolute upper limit on 
compaction parameter, u(= M/a, mass to size ratio of entire configuration in ge- 
ometrized units) < (4/9) for any regular static solution in hydrostatic equilibrium. 2 
Chandrasekhar 3 ' 4 discussed the condition of hydrostatic equilibrium under the 
small adiabatic perturbations, and showed that for each value of the compaction 
parameter, corresponding to the compressible homogeneous spheres, there exists 
a critical (minimum) value of the "ratio of specific heats", j{= Jcrit) such that 
for 7 < jcrit, the configuration becomes dynamically unstable. For the limiting 
case of the compaction parameter approaching the Schwarzschild limit (it = 4/9), 7 
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becomes infinity. 

For dynamically stable superdense objects like neutron stars one may expect a 
finite value of 7. However, for such objects the equations of state are not well known 
[empirically] beyond the density = 10 14 g cm' 3 , 5 and one can only extrapolate the 
equations of state (available in the literature) 6 far beyond this density. As a way 
out, one can impose some restrictions upon the known physical quantities, such 
that, the speed of sound inside the configuration, 



(Where P is the pressure, E is the energy density and s stands for specific entropy) 
does not exceed the speed of light in vacuum, i.e., v < c = 1 (in geometrized units), 
and obtain an upper bound on stable neutron star masses. 7-9 

In the present paper, we have obtained an upper bound on compaction param- 
eter (u < 0.34056 corresponding to a surface redshift of 0.771) for the compressible 
homogeneous spheres, 3 ' 4 by imposing constraint on the "ratio of specific heats", 
7[< 2.59457)], compatible with causality (v < 1) and dynamical stability. This 
value of the compaction parameter is an absolute maximum because, for an as- 
signed value of 7, the maximum compactness would correspond to a compressible 
uniform density sphere, and can be used to obtain an upper bound on neutron star 
masses, as well as the minimum rotation period of a 1.442M© neutron star (the 
maximum mass of the neutron star accurately known at present). 10 

2. Equations Governing Radial Pulsations and Limits on 

Compaction Parameter Imposed by Causality and Ratio of the 
Specific Heats 

Chandrasekhar 3 ' 4 discussed the dynamical stability of fluid spheres with respect to 
small radial adiabatic oscillations on the basis of Einstein's field equations for a 
metric of the form 



where v and A are functions of r and t. By considering the general time-dependent 
field equations appropriate for the metric given by Eq. (1), and letting the quantity 
£(r) represent the amplitude of the 'Lagrangian displacement' from the equilib- 
rium position, namely, £(r, t) = ^{r)e~ lat , where a is the angular frequency of the 
pulsation, the variational base for determining a 2 is given by the equation 3 ' 4 



v = ^(dP/dE) s 



e v dt 2 - e x dr 2 - r 2 d9 2 - r 2 sm 2 0d<f> 2 , 



(1) 
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Causal configurations of homogeneous energy density 3 

A sufficient condition for the dynamical stability of a mass is that the right- 
hand side of Eq. (2) vanishes for some chosen "trial function" £ which satisfies the 
boundary conditions 

£ = at r = 0, (3) 

and 

5P = -£P' - 7 Pe ,y/2 [(r 2 e- ,y/2 £)7r 2 ] 

= at r = a, (4) 

where a is the size of the configuration, 5P is the 'Lagrangian displacement in 
pressure' and the prime denotes radial derivative. The quantity 7 is defined as 3 ' 4 

j=[(P + E)/P](AP/AE) (5) 

where AP and AE denote the 'Eulerian change' in pressure and energy density, 
respectively. With this definition of 7, the boundary conditions [Eqs. (3) and (4)] 
become 

£ = at r = 0, (6) 

and 

AP = -~fPe^ 2 [(r 2 e-^ 2 0'/r 2 } = at r = a. (7) 
For an adiabatic perturbation, using the relation 11-13 

|| = (dP/8E) 8 , (8) 

where s is the specific entropy, and if n denotes the number density, such that, 
P = P(E,n), Eq. (5) becomes 

7 = (n/P)(dP/dn) s = \{P + E)/P}(dP/8E) s , (9) 

or, 

v 2 = *fP/(P + E) = Finite (as long as 7 is finite). (10) 

Let us consider the homogeneous sphere of uniform energy density, E. The equations 
governing equilibrium 3 ' 4 ' 14 can be written in the form of compaction parameter u 
and the radial co-ordinate measured in units of configuration size y(= r/a) as 

87rPa 2 = 6m, 

a n 2 , {l-tuy 2 ) 1 / 2 -(l-Zu) 1 / 2 
8TrPa 2 = 6u v y ' 



3(1 - 2m) 1 / 2 -(1- 2m/ 2 ) 1 / 2 ' 
e,- x = (1 - 2my 2 ), 
e v = (1/4)[3(1 - 2m) 1 / 2 - (1 - 2m/ 2 ) 1 / 2 ] 2 . (11) 
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Table 1. Various values of the compaction pa- 
rameter u and the corresponding critical (mini- 
mum) values of ~f(*f C rit) compatible with dynam- 
ical stability of compressible homogeneous sphere 
of uniform energy density, E. (P/E) and i>o rep- 
resent, respectively, the pressure-density ratio and 
the speed of sound at the center of the configura- 
tion. It is seen that the dynamical stability of causal 
configurations (i. c., u < 0.340555) is constrained 
by the 'ratio of the specific heats', 7 < 2.594570. 



u 


Tcrii 


(P/E)o 


-o 2 


0.050000 


1.38400 


0.02780 


0.03743 


0.100000 


1.44910 


0.06272 


0.08552 


0.150000 


1.53550 


0.10817 


0.14990 


0.200000 


1.65620 


0.17027 


0.24097 


0.250000 


1.83750 


0.26120 


0.38056 


0.277800 


1.98430 


0.33340 


0.49615 


0.300000 


2.14110 


0.40958 


0.62210 


0.320000 


2.32890 


0.50000 


0.77630 


0.333300 


2.49000 


0.57710 


0.91118 


0.340555 


2.59457 


0.62710 


1.00000 



Equation (2) is evaluated for Eqs. (11) with respect to the trial function, 
£ = re u / 4 [because it gives the most rigorous results 15 ' 16 among the various 
trial functions of the form £ = re"/ N ,N = 2, 3, 4. ..00, and the form given as 
£ = 6ir(l + air 2 + a 2 r 4 + a 3 r 6 + ...)e"/ 2 , where ai,a 2 ,a 3 , ... are adjustable 
constants 17 ' 27 ] for various assigned values of the constant 7, so that the config- 
uration becomes compressible, and the speed of sound in this medium remains 
finite and is given by [from Eqs. (10) and (11)] 

v 2 = 7[1 - exp(v - v a ], (12) 

where the subscript 'o' denotes the value of the corresponding quantity at the 
surface of the configuration. 

The compaction parameters of dynamically stable configurations compatible 
with causality [i.e., v < 1] are given in Table 1. These values of 7 are consistent 
with those obtained by Chandrasekhar. 3 ' 4 It is seen that the compactness of the 
causal configuration, u < 0.34056 is constrained by the 'ratio of the specific heats, 
7 < 2.59457. Notice that for a perfectly incompressible homogeneous fluid sphere 
of uniform energy density, E — n, and the ratio of the specific heats, 7, and the 
speed of sound in this medium become infinity for all values of u < (4/9), and the 
configuration would be dynamically stable for all values of u < (4/9). 

3. Application to Obtain an Absolute Upper Bound on Mass and 
Uniform Rotation of Relativistic Stars 

By assigning the energy density E equal to the average nuclear density, i.e. E = 
2 x 10 14 g cm" 3 , 7 in Eqs. (11) for the compaction parameter, u — 0.34056, we obtain 
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Causal configurations of homogeneous energy density 5 

a maximum mass of 5.4M Q [which is larger than the masses obtained earlier]. 7-9 
Although, the assigned value of energy density thus substituted is reasonable, it 
might be fiduciary. Therefore, we have to constrain the value of the energy density 
by some observational fact, and it may be constrained by the fastest rotating pulsar, 
PSR 1937 + 20, with rotation period, P = 1.558 ms, known to date. 18 

The determination of maximum masses of neutron stars corresponding to maxi- 
mum rotation rates require complete general relativistic calculations. 19 ' 20 However, 
the value of E and hence the maximum mass of a static (non-rotating) configura- 
tion, corresponding to higher u values, can be obtained very accurately (see, e. g. 
Refs. 21 and 22) by using the empirical formula given by Koranda, Stergioulas, and 
Friedman 23 in the following form 

JVmm(ms) = 0.740[M ma ,/M o r 1/2 [a ma:c /10km] 3/2 , (13) 

where P ro t,min is the (minimum) rotation period corresponding to a configuration, 
rotating (uniformly) with maximum angular velocity, and M max and a max repre- 
sent, respectively, the maximum mass and the corresponding size of the non- rotating 
configuration. Rewriting Eq. (13) in terms of compaction parameter, u[= M/a], and 
angular velocity Q max (= 2ir / P rot . min ) we obtain 

Q max = 2.21 x \Q w [ulllJa max {cm)]s- 1 (14) 

where u max is the maximum u value of the non-rotating configuration, such that 
the configuration becomes dynamically unstable when u exceeds u max , and a max 
represents the corresponding radius of the configuration. Note that the correspond- 
ing formula which gives an error of 4 - 5 % was previously obtained by Haensel and 
Zdunik. 24 By the use of Eqs. (11) into Eq. (14), we obtain 

£(gcm- 3 ) = 6.59 x lO^JWCs" 1 )] 2 , (15) 

and, 

M max (cm) = [3ul ax /4TTE(cm- 2 )}^ 2 . (16) 

Thus, the uniform energy density of the configuration depends only upon the 
rotation period, and not upon the compaction parameter u. Therefore, it is clear 
from Eqs. (15) and (16) that for a given value of the rotation period P rot the max- 
imum mass of the stable configuration depends only upon the maximum value of 
u(u max ). For P rot — 1.558 ms, Eq. (15) gives the energy density, E, of the config- 
uration as 1.072 x 10 14 g cm~ 3 , the substitution of u max = 0.34056 from Table 
1 into Eq. (16) gives the maximum mass of the configuration, M max = 7.387M Q 
and the corresponding radius, a max = 31.974 km. Notice that the range of energy 
density obtained in this manner is applicable to the baryonic equation of state, 
known as Q-Star equation of state. 5 ' 22 ' 25 

For the value of u = 0.34056, the speed of sound is maximum, v = 1, at the 
center and decreases monotonically from center to the surface of the configuration, 
and at the surface it vanishes along with the pressure. These masses have important 
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implications regarding massive compact objects like Cyg X-l, Cyg XR-1, and LMC 
X-3 etc. These limits of maximum masses could be updated if the faster rotating 
pulsars are observed in the future [see, Ref. 22 for a detailed discussion in this 
regard] . 

On the other hand, if we impose the constraint on energy density [by using Eq. 
(16)] such that the maximum mass of neutron star is 1.442M Q , 10 then Eq. (15) 
gives the minimum rotation period 0.3041 ms for the maximum u = 0.34056, and 
the corresponding uniform energy density E is obtained as 2.813 x 10 15 g cm~ 3 . 

4. Results and Conclusions 

An absolute upper bound on compaction parameter, u < 0.34056 [or the surface 
redshift < 0.771], compatible with causality and the ratio of the specific heats, 

7 < 2.59457, is obtained by using the dynamically stable compressible homogeneous 
sphere. This upper limit of compaction parameter gives 

(i) The maximum static mass of conventional model of neutron stars [taking 
E = 2 x 10 14 g cm~ 3 ] 7 ' 8 as 5.4M . This is greater than 4.8M Q considered as an 
upper limit earlier. 

(ii) The maximum mass of static neutron star exceeds to the value of 7.387M Q 
which is greater than the upper limit of 5.3M Q for neutron stars (so called Q-star 
models) obtained by Hochron, Lynn and Selipesky. 26 This may have important 
implications for the heavy compact objects like Cyg X-l, Cyg XR-1, and LMC-X3 
which may not, necessarily, be black holes. 

(iii) The minimum rotation periods for a static 1.442M neutron star to be 

0. 3041 ms with a uniform energy density E as 2.813 x 10 15 g cm~ 3 . 

(iv) A causally consistent and dynamically stable model of ultra-compact objects 
[u > (1/3)] which has important astrophysical significance [see, e.g. Ref. 16 and 
references therein]. 
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